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Using state-of-the-art first-principles calculations we study the properties of the ferromagnetic 
Heusler compounds NiYSi where Y stands for V, Cr or Mn. NiCrSi and NiMnSi contrary to NiVSi 
are half-metallic at their equilibrium lattice constant exhibiting integer values of the total spin 
magnetic moment and thus we concentrate on these two alloys. The minority-spin gap has the same 
characteristics as for the well-known NiMnSb alloy being around ~1 eV. Upon tetragonalization the 
gap is present in the density of states even for expansion or contraction of the out-of-plane lattice 
parameter by 5%. The Cr-Cr and Mn-Mn interactions make ferromagnetism extremely stable and 
the Curie temperature exceeds 1000 K for NiMnSi. Surface and interfaces with GaP, ZnS and Si 
semiconductors are not half-metallic but in the case of NiCrSi the Ni-based contacts present spin- 
polarization at the Fermi level over 90%. Finally, we show that there are two cases of defects and 
atomic-swaps. The first-ones which involve the Cr(Mn) and Si atoms induce states at the edges of 
the gap which persists for a moderate-concentration of defects. Defects involving Ni atoms induce 
states localized within the gap completely destroying the half-metallicity. Based on single-impurity 
calculations we associate these states to the symmetry of the crystal. 



PACS numbers: 75.47.Np, 71.20.Bc, 71.20.Lp 



I. INTRODUCTION 

Magnetoelectronics, also known as spintronics, is prob- 
ably the field with the most rapid growth in materi- 
als science bringing to the center of scientific research 
new phenomenaii'^'^ For researchers dealing with ab- 
initio calculations one of the most interesting concept 
in spintronics is the half-metallicity.—'^'^ Research has 
been focused initially on the explanation of the origin 
of the half-metallicity and afterwards on the prediction 
of new half-metallic materials with potential applications 
in spintronics devices. Half- metals are hybrids between 
normal metals and semiconductors: the majority-spin 
band is crossed by the Fermi level as in a normal metal 
while the Fermi level falls within a gap in the minority- 
spin band as in semiconductors leading to a perfect 100% 
spin-polarization at the Fermi level^ maximizing the ef- 
ficiency of magnetoelectronic devices<i de Groot and his 
collaborators in 1983 were the first to predict the ex- 
istence of half-metallicity in the case of the intermetal- 
lic semi-Heusler alloy NiMnSbf^ and the origin of the 
gap seems to be well understoodi^ '^^i^^'^^'^'^i^'^ There ex- 
ist several ab-initio calculations on NiMnSb reproducing 
the initial results of de Groot and collaborators , ^^5,16^ 17 
and Galanakis et al. showed that the gap arises from 
the hybridization between the d orbitals of the Ni and 
Mn atoms.— Moreover in the latter reference it was 
shown that the total spin moment for the Heusler al- 
loys of the NiMnSb type follows a Slater-Pauling be- 



havior being in /is the total number of valence elec- 
trons in the unit cell minus 18 since there are exactly 
nine occupied minority-spin statesi^ Theoretical calcula- 
tions for the interfaces of these materials with the semi- 
conductors are few and all results agree that in general 
the half-metallicity is lost both at the surfacesi2iii^'22i'^ 
and the interfaces with binary semiconductorsi^'^'2ii2^ 
Wijs and de Groot have argued than in the case of the 
NiMnSb/CdS (111) contacts the Sb/S interface keeps 
half-metallicity when the S atoms sit exactly on top of 
Sbi^ Moreover taking into account also the reconstruc- 
tion at the surfaces and interfaces can have an impor- 
tant effect on their propertiesi^ Finally we should note 
that several other aspects of NiMnSb have been stud- 
ied using first-principles calculations like the exchange 
constants and Curie temperatures )^'^'2^ the quasiparti- 
cle states^ and the dynamical effects,— the defects and 
doping^'^^i^ the structural stability^i^ the effect of 
spin-orbit coupling^ the fit of model Hamiltonians to 
describe NiMnSb^^'^S orbital magnetismji^ the pressure 
and thermal expansion effects;^ the temperature effeclj^ 
and the magneto-optical properties 4^'^ 

The half-metallic character of NiMnSb in sin- 
gle crystals has been well-established experimen- 
tally. Infrared absorption^ and spin-polarized positron- 
annihilation^i^ gave a spin-polarization of ~100% at 
the Fermi level. High quality films of NiMnSb have been 
grown )^^i^^i^°i^-'^i^^i^'^ but they were found not to repro- 
duce the half-metallic character of the bulk. Values of 
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58% and 50% for the spin-polarization at the Fermi level 
were obtained by Soulen et alM- and by MancofF et al.r^ 
respectively, and recently Zhu et al.— found a value of 
40% using spin-resolved photoemission measurements on 
polycrystalline films. Ristoiu et a/| 57, 58^59,60 glowed that 
during the growth of the NiMnSb thin films, Sb atoms 
segregate to the surface decreasing the obtained spin- 
polarization; they measured a value of ^^30% at 200K, 
while at room temperature the net polarization was prac- 
tically zero. But when they removed the excess of Sb 
by a flash annealing, they managed to get a nearly sto- 
ichiometric ordered alloy surface terminating in MnSb. 
Inverse photoemission experiments at room temperature 
revealed that the latter surface shows a spin-polarization 
of about 67±9% which is significantly higher than all pre- 
vious values There is also experimental evidence that 
for a temperature of '^80 K there is transition from a 
half metal to a normal ferromagnet^^^ but these ex- 
periments are not yet conclusive. Finally, the effect of 
stress on the magnetic anisotropy of thin NiMnSb films 
and the role of defects have been exploredi^i^ 



Based on the success of first-principles electronic struc- 
ture calculations to describe the properties of NiMnSb, 
several authors have predicted new half-metallic Heusler 
alloys crystallizing in the Clb structure of semi-Heusler 
compounds hke NiCrM and NiVM (M= P, As, Sb, S, Se 
and Te)^i^iS and XCrAl (X= Fe, Co, Ni) and NiCrZ 
(Z= Al, Ga, In)i^ Recently, Katayama-Yoshida and col- 
laborators published a paper including also ab-initio cal- 
culations on NiMnSi semi-Heusler alloy, which was pre- 
dicted to have a Curie temperature of 1050 K,— exceed- 
ing even the 730 K shown by NiMnSbi^S This finding 
motivated us to study the electronic and magnetic prop- 
erties of this compound in detail since except the very 
high Curie temperature, it should be easily grown due 
to the well-known and well-controlled growth of NiMnSb 
and it should be compatible with Si (Si crystallizes in 
a diamond structure and the unit cell of NiMnSi is two 
times the unit cell of Si). We decided to expand our study 
to cover also the closely-related NiCrSi and NiVSi alloys. 
NiVSi was found not to be half-metallic at its equilibrium 
lattice constant and thus we focused our study on NiCrSi 
and NiMnSi alloys. In section |TT] we present the details 
of our calculations and in section IIIII the electronic and 
magnetic properties of the bulk phase of these alloys. In 
section IIVI we discuss the robustness of half-mctallicity 
and ferromagnetism, and in section |V] the properties of 
(100) surfaces and of their interfaces with the Si, GaP and 
ZnS semiconductors which have an experimental lattice 
constant almost identical to the theoretical equilibrium 
lattice constants of NiCrSi and NiMnSi. Finally in sec- 
tion |Vl] we discuss the effect of defects and disorder and 
in section fVIII we summarize and conclude. 
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FIG. 1: Schematic representation of the Clj, structure 
adopted by the semi-(or half-)Heusler compounds like NiMnSi 
and of the L2\ structure adopted by the full-Heusler com- 
pounds like Co2MnSi. The lattice consists of 4 fee sublattices. 
The unit cell in the case of semi-Heuslers XYZ is that of a fee 
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lattice with four atoms as basis: X at (0 0), Y at ( 
vacant site at (i i ^) and Z at (| 1 1) in Wyckoff coordinates. 
In the case of full-Heusler compounds the vacant site is also 
occupied by a X atom. 
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FIG. 2: Calculated total energy within the FPLO method vs 
the lattice constant. All three compounds have a theoretical 
equilibrium lattice constant close to 0.54 nm. Note that we 
have scaled for each compound the total energy so that the 
zero energy corresponds to the one of the equilibrium lattice 
constant. 



II. DESCRIPTION OF PRESENT 
CALCULATIONS 

To study the NiYSi alloys we have employed both 
the full-potential nonorthogonal local-orbital minimum- 
basis band structure scheme (FPLO) ''^i^^ and the full- 
potential screened Korringa-Kohn-Rostoker (FSKKR) 
methodf22ii2i Both methods as we will discuss later give 
similar results for the bulk alloys. FPLO is suitable 
to study the disorder within the coherent-potential ap- 
proximation (CPA). FSKKR is suitable to study sur- 
faces and interfaces since its screened character and the 
use of decimation method for the inversion of the Green 
function matrix lead to an order N scaling of the CPU 
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FIG. 3: (Color online)FPLO calculated total density of states 
(DOS) for the NiVSi compound for four different values of 
the lattice constant. Positive values of the DOS correspond 
to the majority (spin-up) electrons and negative values to the 
minority (spin-down) electrons. The zero of the energy has 
been chosen as the Fermi energy. 



time, where N is the number of inequivalent atoms in 
the unit cell. Moreover Dyson equation makes possible 
the study of single impurities within the framework of 
FSKKR in real spaccj ''^'^^ In both cases we employed 
the local-spin-dcnsity approximation (LSDA)21 for the 
exchange-correlation energy within the framework of the 
density functional theoryji^i^^ 

Semi-Heusler alloys (also known as half-Heusler com- 
pounds) crystallize in the Clb structure which consists of 
four fee sublattices and have the chemical formula XYZ 
where X a high-valent transition metal atom. Y a low- 
valent transition metal atom and Z a sp atom. In Fig. 
[T] we present both the Clf, structure of scmi-Hcuslcrs 
and the L2i structure of the fuU-Heusler compounds like 
Co2MnSi. The unit cell is that of a fee lattice with four 
atoms as basis: X atoms occupy the (0 0) site, Y atoms 
the iWj), the (5^5) site is vacant for the semi-Heuslers 
and occupied by X in the case of fuU-Heuslers, and finally 
Z atoms occupy the (| | |) site (we use Wyckoff coor- 
dinates to denote the sites). The zinc-blende structure 
adopted by a large number of semiconductors, like GaP 
and ZnS is also consisting of four fee sublattices. In the 
case of GaP the (000) site is occupied by a Ga atom, the 
(i i i) site by a P atom, while the other two sites are 
unoccupied. The diamond structure of Si is also coherent 
with the Clb structure of semi-Heusler alloys since dou- 
ble the unit cell gives the zinc-blende structure of binary 
semiconductors. This close structure similarity makes 
the Heusler alloys compatible with the existing scmicon- 



TABLE I: (Color online) Atom-resolved and total spin mag- 
netic moments in /is calculated using both the FPLO and 
FSKKR methods. In the last two columns we present our 
results on the Curie temperature in Kelvin units calculated 
with both the mean-field {T^'^^) and random-phase (Tq 
approximations. 
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FIG. 4: (Color online)Ni and Cr(Mn) resolved DOS for both 
NiCrSi and NiMnSi compounds using both the FPLO and 
FSKKR methods. In the insets we present the DOS for a 
wider energy range. 



ductor technology and thus very attractive for industrial 
applications. 

Since the compounds under study have not been grown 
experimentally, we have to determine their equilibrium 
lattice constant. In Fig. [2] we plot the FPLO calculated 
total energies vs the lattice constant for all three NiVSi, 
NiCrSi and NiMnSi. We have assumed that all three 
arc ferromagnets (which is confirmed by our calculated 
exchange constants, at least for NiCrSi and NiMnSi, in 
section IIV[) . We have scaled the total energy for each 
compound in a way that the zero of the total energy 
corresponds to the minimum energy at the theoretical 
equilibrium lattice constants. All three compounds have 
an equilibrium lattice constant very close to 5.4 A. We 
have decided to use the latter one for all calculations 
presented in this manuscript. 

We should also discuss here the reason that led us to 
neglect NiVSi from our study. In Fig. [3] we present the 
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total DOS for NiVSi and for four different lattice con- 
stants: 5.2, 5.4, 5.6 and 5.8 nm (we remind here that 1 
nm — 10 A). There is a gap both in the majority and 
minority spin bands and the Fermi level for the equi- 
librium lattice constant crosses the conduction minority 
spin-band and the alloy is not half- metallic. Moreover 
it is located in a peak of the majority-spin DOS mak- 
ing fcrromagnctism unstable and probably the compound 
prefers to be non-magnetic. The calculated total spin 
moment presented in Table U is 0.874 /i^ lower than the 
ideal 1.0 hb requested by the Slater-Pauling behavior 
for half-metallic semi-Heusler alloysi^ We have not cal- 
culated the difference between the non-magnetic and fer- 
romagnetic states since it exceeds the scope of this paper. 
Although the Fermi shifts lower in energy as we expand 
the lattice and the compound becomes half-metallic the 
extremely high majority DOS at the Fermi level persists 
and ferromagnetism seems to be unstable. The latter 
argument is also confirmed by the calculated exchange 
constants presented in Fig. [71 The Ni-V interaction is 
negligible and the V-V exchange constant between vana- 
dium nearest-neighbors presented in the figure is very 
weak being 8 times smaller than the Mn-Mn interaction 
in NiMnSi and thus the tendency to ferromagnetism is 
not very strong. Also the Curie temperature deduced us- 
ing these exchange constants is ~200 K within both the 
mean-field and random-phase approximations and thus 
it is considerably smaller than the room temperature at 
which realistic devices should operate. Since NiVSi is not 
so interesting for applications with respect to NiCrSi and 
NiMnSi, we have decided not to include it in the rest of 
our study. 



III. ELECTRONIC AND MAGNETIC 
PROPERTIES 
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FIG. 5: (Color online) FSKKR-calculated spin-resolved band 
structure for the NiMnSi compound along several high- 
symmetry axis. We have denoted the character of the bands 
at the r point and in the brackets we denote the orbitals 
transforming following each representation. 



We will start our discussion from the electronic and 
magnetic properties of NiCrSi and NiMnSi establishing 
also the equivalence between the FPLO and FSKKR 
methods. In Fig. |3]we present the DOS of the transition- 
metal atoms Ni and Cr(Mn) since the DOS of Si is negli- 
gible within both FPLO (dashed hue) and FSKKR (solid 
line). Both methods give almost the same DOS and the 
only difference is the intensity of some peaks. Especially 
concerning the gap (in the right panel we have blown- 
up the region around the Fermi level) both FPLO and 
FSKKR give the same gap-width, the same position of 
the Fermi level with respect to the edges of the gap and 
the same majority-spin DOS at the Fermi level. In Table 
Uwc present also the spin moments using both electronic 
structure methods. FPLO gives integer values of the to- 
tal spin moments while FSKKR gives slightly smaller val- 
ues due to the i summation used in FSKKR. This prob- 
lem has been extensively discussed in Ref. H^. FSKKR 
overestimates with respect to FPLO the hybridization 
between the Ni and Mn d-states. Since FPLO gives a 
more atomic-like description of the Cr and Mn, it pre- 



dicts larger Cr and Mn spin moments and consequently 
smaller Ni spin moments to keep constant the total spin 
moment. Si spin moment shows also a small variation 
between the two methods since its moment arises from 
the hybridization with the Ni-Cr(Mn) d hybrids. The 
total spin moment does not exactly equals the sum of 
the atomic spin moments of the three constituent atoms 
since also the vacant site carries a very small spin mo- 
ment which we do not present in the table. 

Since we have established the cquilivalcnce between 
the FPLO and FSKKR we will go on with the discus- 
sion of the electronic and magnetic properties of NiCrSi 
and NiMnSi. NiMnSi has 21 valence electrons per unit 
cell, one less than NiMnSb, and NiCrSi 20 valence elec- 
trons. In accordance with the Slater-Pauling rule for the 
perfect half-metallic semi-Heusler compounds presented 
in Ref. the calculated total spin moment in the unit 
cell is for NiMnSi 21-18= 3 fiB and 20-18= 2 for 
NiCrSi. The spin moment is mainly carried by the low- 
valent transition metal atoms Cr and Mn. Ni carries a 
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TABLE II: (Color online)FPLO calculated width of the 
minority-spin gap in eV for the cases presented in Fig. [6]with 
respect to the perfect cubic lattice (the ration equals to zero 
since c = a). Negative values of the ratio in the first column 
correspond to contraction and positive values to expansion. 





Gap- Width (eV) 


c — a 
a 


NiCrSi 


NiMnSi 


-5% 


1.029 


0.966 


-2% 


1.067 


1.050 


-1% 


1.071 


1.055 


0% 


1.071 


1.055 


1% 


1.021 


1.017 


2% 


0.983 


0.987 


5% 


0.882 


0.903 



small spin moment of about 0.1-0.2 due to the hy- 
bridization with the Cr(Mn) d-states. Si shows a small 
spin moment which is antiparallel to the moment of the 
transition-metal atoms since the majority p-statcs of Si 
arc more extended in energy and thus partially unoccu- 
pied leading to a surplus of minority-spin Si p-states and 
a negative atomic spin moment. 

The origin of the gap is similar to the case of NiMnSb.— 
Ni and Cr(Mn) d-states hybridize creating 5 bonding d 
hybrids and 5 antibonding ones per spin-direction. Si 
provides one s- and three p-states per spin low in en- 
ergy. The p-states accommodate also d electrons of the 
transition metal atoms. In Fig. [5] we present the band- 
structure along several high symmetry directions for both 
spin-directions for NiMnSi. The character of each band 
at the r point reveals the character of the associated or- 
bitals in the real space. First for the minority-spin band 
structure there is one s and a triple-degenerated at the F 
point p-band coming from the Si atom. Then follow the 
double-degenerated Cg and triple-degenerated t2g bands 
due to the Ni-Mn bonding d hybrids. The t2g band is 
separated by an indirect V ~ X gap from the antibond- 
ing bands. Thus in total there are exactly nine occupied 
minority-spin bands. In the majority spin band we have 
12 instead of 9 atoms. The two out of the three extra 
electrons occupy the antibonding Cg states and the last 
electron occupies partially the antibonding t2g bands. In 
the case of NiCrSi the majority-spin band has to accom- 
modate 11 electrons and not 12 and thus it is enough 
to occupy the majority-spin antibonding eg states. The 
band associated with the latter states is very narrow and 
is also well separated in energy from the antibonding t2g 
bands and this is why in the case of NiCrSi the Fermi level 
is located in a very narrow deep separating two peaks in 
the Cr majority-spin DOS. 

IV. ROBUSTNESS OF HALF-METALLICITY 
AND FERROMAGNETISM 

We will continue our study discussing the robustness 
of the half-metalicity and ferromagnetism. Since alloys 
grow on top of semiconductors, the lattice mismatch can 
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FIG. 6: (Color online)FPLO total DOS of NiCrSi and 
NiMnSi compounds under 1% and 5% expansion and con- 
traction of the c lattice parameter along the 2-axis resulting 
to a tetragonalization of the lattice. The in-plane lattice pa- 
rameter was kept constant and equal to 0.54 nm. 



lead to a tetragonalization of the lattice of the Heusler 
alloys to minimize the elastic energy of the film. Thus we 
have to study for realistic applications the effect of tetrag- 
onalization on the half-metallicity. In Fig. [5] we present 
the total DOS for NiCrSi and NiMnSi using the FFLO 
code under contraction or expansion of the out-of-plane 
lattice parameter by 1 and 5% keeping the in-plane lattice 
constants equal to 0.54 nm. The latter value of tetrago- 
nalization is very large and in most realistic applications 
the lattice relaxes much less. As can been deduced from 
the figure contraction leads to a marginal shift of the gap 
while expansion leads to a very small narrowing of the 
gap. To give more details on the behavior of the gap we 
present in Table [ll] the width of the gap for the different 
cases. For the perfect cubic NiCrSi and NiMnSi alloys 
the gap is 1.071 and 1.055 eV, respectively. Contraction 
by 5% leads to gap-width of 1.029 and 0.966 eV respec- 
tively, while expansion by 5% to 0.882 and 0.903 eV for 
NiCrSi and NiMnSi respectively. Thus tetragonalization 
leaves the gap almost unchanged and since for both al- 
loys the Fermi level is near the middle of the gap this has 
no effect on the half-metallic properties of the two alloys. 

Stable half-metallicity does not mean stable ferromag- 
netism. To study how robust is the ferromagnetism we 
have to calculate the exchange constants and then use 
them to calculate the Curie temperature. We have used 
the formalism presented in Ref. [23 and in Fig. [7] we 
present the Cr-Cr and Mn-Mn exchange constants as a 
function of the distance between the atoms. We neglect 
the Ni-Cr and Ni-Mn interactions since they are very 
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FIG. 7: (Color online) Calculated exchange constants between 
the low-valent transition metal atoms sitting at the Y site (V, 
Cr or Mn) as a function of the distance between them. 



small even if each Cr or Mn atom has four Ni atoms 
as first neighbors. The much larger Cr and Mn spin mo- 
ments are responsible that interactions between them are 
much larger than between Ni and Cr(Mn) atoms. In the 
case of NiMnSi the exchange constants are very large 
and even Mn-Mn next-nearest neighbors contribute to 
the stability of ferromagnetism. The exchange constant 
between nearest neighbors, usually in literature denoted 
as Ji, is 1.2 niRy. This is a very large value and leads to 
very large values for the Curie temperature (Tq ). Cal- 
culating the latter one as explained in Ref. |23 we get 
a value of 1505 K within the mean field approximation 
(MFA) and 1129 K within the random-phase approxima- 
tion (RPA) (see Table [I| close to the calculated value of 
1050 K in Ref. [6^. Our experience on NiMnSb shows 
that the RPA results for semi-Heusler compounds are 
closer to the experimental values than MFA. In the case 
of NiCrSi only the interaction between Cr atoms nearest- 
neighbors contributes to the stability of ferromagnetism 
and although its strength is smaller than the Mn-Mn one, 
it is enough to get a Curie temperature as high as 900 K 
within MFA and 698 K within RPA. Thus both NiCrSi 
and NiMnSi are very promising for realistic applications 
since their Tq exceed by far the room temperature of 
~300 K. 



V. SURFACES AND INTERFACES 
A. (001) Surfaces 

Although surfaces are not as interesting as interfaces 
we will start our discussion from the (001) surfaces to 
give a basis for our discussion in the case of interfaces. 
We have chosen this surface orientation since this is the 
most usual one encountered in experiments. We have to 
note first that we have considered only abrupt surfaces 
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FIG. 8: (Color online)FSKKR-calculated atom-resolved DOS 
for the Ni and Cr(Mn) atoms sitting at the first two surface 
layers for the NiCrSi(OOl) (left panel) and NiMnSi(OOl) (right 
panel) surfaces. Note that e.g. there are two different termi- 
nations for the NiCrSi(OOl) surface: a Ni surface layer with a 
CrSi subsurface layer or a CrSi surface layer with a Ni sub- 
surface layer. 



and have not taken into account possible reconstructions. 
When we open the (001) surfaces in NiCrSi and NiMnSi 
we have two possible terminations. The surface layer 
can be either a Ni-Void one with a Cr(Mn)Si as subsur- 
face layer, or vice versa. For the completeness of our 
results we have considered both terminations and have 
used a slab of 15 layers embedded in vacuum which has 
been shown to be adequate for the half-metallic Heusler 
alloys (see Ref. [l^. We have also to note here that 
the Ni-terminated surfaces contain also a vacant site and 
are susceptible to show large reconstructions with re- 
spect to the CrSi- and MnSi-terminated NiCrSi(OOl) and 
NiMnSi(OOl) surfaces, respectively]^ Thus we will con- 
centrate mainly on the latter surfaces. 

In Fig. [5] we present the DOS of the transition metal 
atoms at the surface and subsurface layers for both 
NiCrSi and NiMnSi alloys. It is obvious that in all cases 
the surfaces are not half- metallic. The bonds which are 
broken at the surface change the energy position of the 
states and the gap vanishes similarly to what happens 
to semiconductors where surface states make their sur- 
faces metallic in most cases. More information we can 
get from Table IIIII where we present the spin moments 
of the atoms at the surface and subsurface layers (which 
are significantly enhanced with respect to the bulk due 
to the missing hybrization) and also the number of ma- 
jority (N^) and minority (N^) electronic states at the 
Fermi level for these atoms and the spin-polarization, P, 
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TABLE III: (Color online)Upper panel: FSKKR calculated atomic spin moments for the NiCrSi(OOl) and NiMnSi(OOl) surfaces 
for both possible terminations. We present the moments for the atoms in the first two surface layers and in the last column 
their sum. Lower panel: for each one of the atoms in the upper panel we present the ratio of the majority-spin DOS at the 
Fermi level (A'^^) vs the minority-spin DOS at the Fermi level [N^). In the last column we present the spin-polarization at the 
Fermi level calculated as .^"'Y taking into account both the surface and subsurface layers, P. 
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defined as taking into account both the surface 

and subsurface layers. 

We will start our discussion from the CrSi- and MnSi- 
terminated surfaces. In the case of the CrSi-terminated 
surfaces the DOS with the exception of the gap area is 
very similar to the bulk calculations. The Ni atom in the 
subsurface layer presents practically a half-ferromagnetic 
character with an almost zero spin-down DOS, while for 
the bulk NiCrSi there is an absolute gap. The Cr and Si 
atoms in the surface layer show more pronounced differ- 
ences with respect to the bulk, and within the gap there 
is a very small Cr-d DOS. These states are strongly lo- 
calized at the surface layer as at the subsurface layer 
there is practically no states inside the gap. When we 
open the CrSi(OOl) surface, the Cr atoms at the sur- 
face layer loose two out of their four nearest Ni neigh- 
bors. This breaking of symmetry at the surface pushes 
the Cr majority-spin states deeper in energy to account 
for the missing hybridization and the charge accommo- 
dated in these bonds is now accommodated by the Cr 
majority-spin states leading to a considerable increase 
of the Cr atomic spin moment as it can be seen also in 
Table IIIII Moreover the splitting between the unoccu- 
pied Cr states above Ep and the center of the occupied 
Cr states decreases and at Ep a surface state appears. 
This behavior of the CrSi-terminated NiCrSi(OOl) sur- 
face is in agreement with previous calculations on the 
MnSb-terminated NiMnSb(OOl) surfaceJ^ Similar is the 
situation also for the MnSi-terminated NiMnSi(OOl) sur- 
face. Mn atoms act like the Cr ones and accommodate 
most of the extra-charge due to the broken-bonds in the 
majority-spin band. But Mn atoms have one valence 
electron more than Cr atoms which in the bulk is al- 
ready accommodated in the majority band and thus for 
Mn atoms the energy cost to accommodate extra charge 



in the majority band is larger than for the Cr atoms. As 
a result part of the extra-charge due to the broken bonds 
is accommodated also in the minority-spin band leading 
to a smaller increase of the spin moment with respect to 
the Cr atoms and a more intense surface state within the 
gap. Mn atoms at the surface layer show an increase of 
their spin moment of about 0.6 with respect to the 
bulk reaching a value of '^3.4 |J,B^ while in the case of Cr 
atoms the increase of the spin moment with respect to 
the bulk is 0.9 ^b- 

In the case of the Ni terminated NiMnSi(OOl) surface, 
the Ni atoms at the surface layer show a largely differ- 
ent shape of the DOS with respect to the Ni atom at 
the subsurface layer in the case of the MnSi-terminated 
surface and the bands are located higher in energy with 
respect to the MnSi(OOl) termination. The Mn atoms 
at the subsurface layer show DOS very similar to the 
bulk case and the spin magnetic moment is ~3.17 ^b 
very close to the bulk value of ~2.84 ^b- In this case it 
is clear that there is a minority surface state located ex- 
actly at the Fermi level since a very narrow peak centered 
at the Fermi energy appears in the minority-spin DOS. 
This surface state survives also for the Mn atoms at the 
subsurface layer. Similar arguments are also valid for the 
Ni-terminated NiCrSi(OOl) surface. Cr atoms behave as 
in the bulk with slightly larger spin moments. The spin 
magnetic moment of the Ni atoms at the surface layer is 
considerable larger than the value for the bulk (0.28/iB 
for the Ni surface atom in NiCrSi with respect to 0.17 
jjLB for the Ni atom in the bulk NiCrSi, and 0.46/iB for 
the Ni surface atom in NiMnSi with respect to 0.23 ^b 
for the Ni atom in the bulk NiMnSi) The surface state in 
the case of the Ni-terminated NiCrSi(OOl) surface is more 
broad than in the case of the NiMnSi alloy and cannot 
be distinguished from the rest of the DOS. 
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For applications more important is the spin- 
polarization at the Fermi level. We present in the lower 
panel of Table IIIII the number of states at the Fermi 
level for each atom and for both spin directions. We 
calculate the spin polarization as the number of major- 
ity DOS divided by the total number of electronic states 
at the Fermi level. In the case of Ni-tcrminated NiMnSi- 
surfaces the number of states for the two spin directions is 
equal due to the surface states and the spin-polarization 
is around 50% while for the MnSi-terminated surface the 
surface state is less intense and the spin polarization 
reaches 67% . In the case of NiCrSi surface Cr atoms 
show a very large majority-spin DOS at the Fermi level 
with respect to the Mn atoms since the Fermi level crosses 
the eg states. This washes out the weight of the surface 
states and NiCrSi surfaces show very large values of the 
spin polarization: 82% for the Ni-terminated surface and 
86% for the CrSi-terminated surface. These results do 
not mean that NiCrSi is more suitable for applications 
since the electronic properties change considerably at the 
interfaces with semiconductors. 



B. Interfaces 

We will continue our study with the interfaces between 
the half-metal and the semiconductors. Such structures 
are assumed to present accurately the realistic situation 
since in most devices half-metals are employed to inject 
current into a semiconductor. We have assumed that the 
stacking direction is the (001) and used a slab with 15 
layers of the half-metal and 9 layers of the semiconductor 
in order to have two equivalent interfaces as in Ref. [23l . 
Moreover we have searched for semiconductors crystalliz- 
ing in the zinc-blende structure which is compatible with 
the lattice structure of the Heusler alloys^S and with a 
lattice constant close to the 5.4A which is the equilibrium 
lattice parameter of both NiCrSi and NiMnSi. Such semi- 
conductors are the binary III-V semiconductor GaP, the 
II- VI semiconductor ZnS and finally the IV semiconduc- 
tor Si. We have considered all possible interfaces and in 
Table IIVI we resume our results by presenting the spin- 
polarization at the Fermi level taking into account the 
first two half-metallic interface layers, Pi, and both the 
first two half-metallic and the first two semiconducting 
layers, P2- Note here that contrary to previous studies 
we define the spin-polarization P as the ratio between 
the number of majority (spin-up) states and the total 
number of states at the Fermi level and thus P repre- 
sent directly the proportion of the majority electrons at 
the Fermi level. There are two general remarks in this 
table. First, when we consider also the semiconductor 
layers the spin-polarization decreases slightly by a few 
percent in all cases. The interface half-metallic layers 
polarize the semiconducting layers at their vicinity and 
the first few semiconducting layers become metallic. But 
the induced DOS of the semiconductor within the gap 
is similar for both spin-directions decreasing the spin- 
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FIG. 9: (Color online)FSKKR-calculated atom-resolved DOS 
for the Ni and Cr atoms at the interface for the following cases: 
(i) upper left panel: Ni/Ga structure for the NiCrSi/GaP 
interface, (ii) upper right panel: Ni/P structure for the 
NiCrSi/GaP interface, (iii) middle left panel: Ni/Zn structure 
for the NiCrSi/ZnS interface, (iv) middle right panel: Ni/S 
structure for the NiCrSi/ZnS interface, (v) lower left panel: 
Ni/Si structure for the NiCrSi/Si interface, and (vi) lower 
right panel: CrSi/Si structure for the NiCrSi/Si interface. 

polarization. Second, the interfaces between NiMnSi and 
the three semiconductors show considerably smaller spin- 
polarization than the corresponding interfaces created by 
NiCrSi. This is expected since, as for the surfaces, the 
Cr atoms have a very high majority DOS at the Fermi 
level with respect to the Mn atoms and they minimize 
the effect of interface states. Thus in the following we 
will concentrate on the interfaces between NiCrSi and 
the three different semiconductors. 

When the contact at the interface is created by the Ni 
atomic layer and not the CrSi layer the obtained spin- 
polarization is much higher and can exceed 90% as in 
the case of Ni/P and Ni/Si contacts. The different be- 
havior between the Ni and CrSi interface layers can be 
easily understood in terms of symmetry and hybridiza- 
tion. Ni interface atoms have a spin moment close to 
their bulk value as shown in Table |V] where we present 
the atomic spin magnetic moments for all atoms at the 
interface. In the bulk case Ni has 4 Cr and 4 Si atoms 
as first neighbors. On the Ni-terminated (001) surface 
the Ni atom loses half of its first neighbors. But if an 
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TABLE IV: (Color online) Similar to the lower panel of Table ITITl for all studied interfaces between the half- metals NiCrSi and 
NiMnSi and the GaP, ZnS and Si semiconductors, taking into account all possible interface structures. The spin-polarization 
ratios at the Fermi level have been calculated as jy/^^^j^ taking into account only the two first interface layers of the half-metal, 
Pi, and the first two interface layers of the half-metal and the first two interface layers of the semiconductor. Pi- 
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TABLE V: (Color online)FSKKR calculated atomic spin mo- 
ments for the two layers of the NiCrSi half-metal and of the 
two first layers of the GaP, ZnS or Si semiconductors at the 
interface. Void denotes the empty site in the same layer with 
the Ni atom. Vcl denotes the empty site in the same layer 
with Ga(Zn) and Vc2 in the same layer with the P(S) atoms. 
In the case of the interfaces with Si, we denote with Sil the 
Si atoms at the interface layer and with Si2 the Si atoms at 
the subinterface layer. 
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interface with a semiconductor is formed, the nickels two 
lost Si neighbors are replaced by two sp atoms presenting 
similar electronic structure to the Si ones, although they 
have in general different number of valence electrons, and 
the situation is much closer to the bulk case than when 
we have CrSi as an interface layer. Now the Si p bands at 
lower energy are not destroyed since the sp atoms have a 
behavior similar to Si and still they accommodate tran- 



sition metal d electrons. Thus the only change in the Ni 
DOS comes from the missing two Cr neighboring atoms. 
The DOS of the Ni atom at the interface layer presented 
in Fig. [9] for several Ni/sp contacts is clearly very close 
to the bulk case. Moreover the Cr atoms at the subin- 
terface layer are also almost bulk like showing with the 
exception of the gap area a DOS close to the bulk one for 
all cases of Ni/sp contacts. Between the different cases 
presented in Fig. [9] -Ni/Ga and Ni/P contacts in the 
case of NiCrSi/GaP interfaces, Ni/Zn and Ni/S contacts 
in the case of NiCrSi/ZnS interfaces and Ni/Si contact 
in the case of NiCrSi/Si interface- the DOS of the Ni in- 
terface atoms and Cr subinterface atoms only marginally 
change. The small variations at the gap area arise from 
the different hybridization between the t2g d-orbitals of 
Ni and the p-orbitals of the sp atom which transform with 
the same representation and thus can couple. When the 
sp atom is Si the situation is closer to the bulk NiCrSi 
and we get a spin polarization as high as 93%. When 
the sp atom is P which has one electron more than Si 
the P2 drops slightly to 90%. As the difference between 
the electronic structure of Si and the sp atoms becomes 
more important the hybridization effects become more 
pronounced inducing a larger DOS in the minority-spin 
band and thus leading to smaller values of P2- 

In the case when the contact is made up by the CrSi 
layer the Cr atoms have an immediate environment very 
different from the bulk NiCrSi and thus their electronic 
properties are considerably altered with respect to the 
bulk case. In Fig. [9] we present the DOS of the Cr atom 
at the interface layer and the Ni atom at the subinterface 
layer in the case of a CrSi/Si contact. The majority-spin 
bands are broader with respect to the Ni/Si contact and 
they are pushed deeper in energy. This attracts some of 
the minority unoccupied Cr states also deeper in energy 
and now the Fermi level crosses the antibonding minor- 
ity d-bands leading to a smaller spin-polarization at the 
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Fermi level. While the CrSi-based contacts keep a degree 
of polarization close to 60% due to the large intensity of 
the Cr majority-spin DOS at the Fermi level, in the case 
of the MnSi-based contacts the situation is even worse 
and P2 becomes smaller than 50% meaning that at the 
Fermi level the spin-down states are more than the spin- 
up states. 

Finally we should also briefly discuss the spin mag- 
netic moments presented in Table |V] for the case of 
NiCrSi/(GaP or ZnS or Si) interfaces. The atoms in 
the semiconductor spacer become magnetic due to the 
polarization from the transition-metal atoms in the half- 
metallic spacer but their spin moments arc very small 
(less than O.OI^b in all cases). Ni and Cr atoms in all 
cases present spin moments close to the bulk values and 
there is no general trend associating the spin magnetic 
moments to the spin-polarization of the interfaces. Thus 
for each case we have to study the DOS at the inter- 
face to get valid conclusions about the behavior of the 
half-metallicity. 

We should also note here that we have neglected re- 
laxation effects at the interface. Relaxation effects mod- 
ify the spin-polarization at the Fermi level since it leads 
to increased hybridization and charge-transfer. It was 
shown in the case of a Ni/P contact in the case of 
the NiMnSb/GaP interface that the Ni-P interlayer dis- 
tance was reduced by 18%, while the neighboring inter- 
layer distances where expanded by 5-7% as compared to 
the ideal bulk values leading to a decrease of the spin- 
polarizationi^ Nevertheless, these effects did not destroy 
the gap but shifted the local DOS slightly deeper in en- 
ergy and a suitable doping at the interface restored the 
spin-polarization. Moreover in the case where the spin- 
polarization for the unrelaxed surface is as high as 90% 
as in some of the cases presented above, the effect of 
relaxation will be much weaker. 



VI. DISORDER AND DEFECTS 

Our final section is devoted to the effect of defects. 
Defects in large concentrations can induce new bands 
within the minority-spin gap and completely destroy half- 
metallicity. But even in low concentrations they can cou- 
ple to interface states and the latter ones become con- 
ducting destroying the spin polarization during transport 
experiments through a half- metal/semiconductor inter- 
face. To perform our study we have employed the FPLO 
method within the CPA approximation to account for 
disorder in a random way. We have studied two types of 
defects (i) the creations of Mn(Cr) and Si antisites and 
Mn(Cr)-Si swaps which keep the half-metallic character, 
and (ii) the creation of Ni antisites and Ni-Mn(Cr) atomic 
swaps which induce states within the gap destroying the 
gap. 



TABLE VI: (Color online)FPLO calculated width of the 
minority-spin gap in eV for the cases presented in figure [TITl as 
a function of the disorder concentration x. With "-" we de- 
note the cases where the gap persists but the half-metallicity 
is lost since the Fermi level is above the gap. 
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A. Defects conserving the half-metallicity 

As we have mentioned just above the creation of 
Cr(Mn) and Si antisites keeps the half- metallic charac- 
ter, as-well-as the Cr(Mn)-Si atomic swaps. The reason 
is pretty simple. If we examine the band-structure of the 
minority-spin electrons in Fig. [5]the states just below the 
gap are the t2g d-hybrids which are mainly located at the 
Mn(Cr) sites. Above the gap we have the antibonding eg 
d-states and then the antibonding t2g d- hybrids. The t2g 
states transform using the same representation as the p 
states of the Si atom. The sites occupied by the Si and 
Cr(Mn) atoms in the perfect compounds have the same 
symmetry, with four Ni and four Void sites as first neigh- 
bors, rotated by 90°. Thus the equivalence between the 
t2g and p states is obvious. When Cr(Mn) atoms mi- 
grate to Si sites or Si atoms to Cr(Mn)sites, it results 
to a broadening of the bonding t2g hybrids just below 
the gap. Also the antibonding t2g band above the gap 
broadens and it overlaps with the antibonding Cg states. 
These two simultaneous phenomena lead to a shrinking 
of the gap which as we will discuss below persists for 
a moderate percentage of disorder. Of course above a 
critical value of disorder the gap vanishes. In the case 
of Cr(Mn)-Si atomic swaps, the phenomenon discussed 
above is more intense and the gap vanishes for a smaller 
degree of disorder. 

In Fig. [To] we present the total, Ni, and Cr(Mn) DOS 
for the case of 10% Cr(Mn) antisites [NiCr(Sio.9Cro.i) 
and NiMn(Sio.9Mno.i) alloys], the case of 10% 
Si antisites [Ni(Cro.9Sio.iSi) and Ni(Mno.9Sio.i)Si], 
and 10% atomic swaps [Ni(Cro.9Sio.i)(Sio.9Cro.i) and 
Ni(Mno.9Sio.i)(Sio.9Mno.i)]. We also present for the first 
and third case the DOS of the Cr(Mn) impurity atoms 
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FIG. 10: (Color online) Left panel: FPLO-calculated total and atom resolved DOS for three distinct doping cases in NiMnSi 
which conserve the half-metallicity. Right panel: similar DOS for NiCrSi. Atomic DOS have been scaled to one atom. 



at the antisitcs. First remark is that Ni DOS is similar 
for all three compomids and similar to the Ni DOS in 
the perfect compounds. The same is true for the Cr and 
Mn atoms at the perfect sites. Thus the antisites and 
atomic swaps have only marginal effect on the properties 
of the transition-metal atoms at perfect sites. Concern- 
ing now the impurity Cr and Mn atoms sitting at perfect 
Si sites, the hybridization between the Cr(Mn) t2g states 
and p states leads to important differences with respect 
to the Cr and Mn atoms sitting at the perfect sites. The 
minority-spin impurity states move lower in energy but in 
the case of the Mn impurities the large exchange splitting 
between the majority and minority d-states keeps these 
states well above the gap. In the case of the Cr impurities 
the exchange splitting is significantly lower than for the 
Mn impurities and the Cr unoccupied minority eg states 
are crossed by the Fermi level in the case of atomic swaps 
(lower panel). In the case of Cr antisites the antibond- 
ing Cg states move even lower in energy and the Fermi 
level falls within a tiny gap separating the antibonding eg 
and t2g states and the Cr impurity atoms have an almost 
vanishing spin magnetic moment. (We do not present 
the spin moments of the other atoms since they are very 
close to the values for the perfect compound.) 

To confirm the character of the Cr impurity states we 
have also employed FSKKR method and have made cal- 
culations for a single Cr impurity in NiCrSi and in Fig. 
[TT] we present the d states of Cr antisites at a Si site 
projected on the eg and t2g orbitals. For the impurity 



calculations we perform the self-consistent calculations 
for NiCrSi and produce the Green function for the bulk 
compoundi ^^i^^ Then we consider a cluster of 65 atoms 
surrounding the Si site embedded in an infinite NiCrSi 
bulk crystal and calculate the electronic structure of the 
atoms in this cluster in real space considering a Cr atom 
at the central Si site. (We have to note here that we 
have neglected the relaxation of the atomic positions 
within the cluster due to the impurity atom, which is 
very tedious computationally, since we expect them to 
be small due to the high symmetry of the crystal and 
an extensive study of impurities is out of the scope of 
the present manuscript.) As we can deduce from the 
DOS in Fig. [5] the antibonding unoccupied minority d- 
statcs have mainly their weight at the Cr site while the 
bonding occupied minority states at the Ni site. The mi- 
nority Cg-states of the Cr-impurity are located lower in 
energy with respect to the Cr atom in the perfect site. 
Moreover they are located exactly above the Fermi level 
and are also very localized in energy. The t2g minority 
states of the Cr-impurity atom form a wider band which 
starts just at the right edge of the eg peak and goes very 
high in energy. In the case of atomic swaps the eg states 
move even lower in energy forming this tiny gap. The Cr- 
impurity atom polarizes slightly the neighboring atoms 
but the effect very quickly vanishes inside the cluster. As 
we can deduce from these results on a single impurity in 
a perfect bulk crystal, which are in very good agreement 
with the results obtained for a disordered crystal within 
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FIG. 11: (Color online) d-resolved atomic DOS for several 
cases of single impurities in NiCrSi compound calculated us- 
ing the FSKKR method. The DOS has been decomposed 
in the double-degenerated eg and triple-degenerated t2g con- 
stituents. 



the coherent potential approximation (CPA) where each 
disordered site is occupied by a pseudoatom with mixed 
Si-Cr properties, local-order effects and long-range effects 
have the same influence on the electronic properties of 
NiCrSi alloy. Wc have performed similar calculations for 
a single Mn impurity atom sitting at a Si site in NiMnSi 
and the results were similar to the case of a Cr antisite in 
NiCrSi but now the Cg and t2g states in the minority-spin 
band where higher in energy due to the larger exchange- 
splitting of the Mn occupied majority and unoccupied 
minority d states. 

Finally, wc should discuss the width of the gap which 
is the most interesting property affected by the defects. 
In tablc lVll we present the width of the gap as a function 
of the concentration of defects for all three studied cases 
and for both Cr- and Mn-based compounds. The more 
mild effect is the creation of Si antisites (middle panel of 
the table) since no-transition-metal impurities exist. The 
gap persists with large width-values even for as much as 
20% of Si antisites (0.626 eV for the Mn-based compound 
and 0.613 cV for the Cr-based alloy). When we create 
atomic swaps (lower panel of the table) the gap-width 
shrinks with a double rate and for NiMnSi alloys it drops 
to 0.216 eV for the case of 20% atomic swaps. In the case 
of NiCrSi the gap-width is half its initial value for 5% of 
atomic swaps (0.484 eV instead of 1.071 eV). For 10% 
creation of Cr-Si atomic swaps the gap persists but the 
Fermi level is above the gap and the defected alloy is not 
half-metallic. As we have mentioned above in the case 
of simple Mn(Cr) antisites the gap shrinks very quickly 
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FIG. 12; (Color online)Left panel: FPLO-calculated to- 
tal and atom resolved DOS for two distinct doping cases in 
NiMnSi which destroy the half-metallicity. Note that with E 
we denote the vacant site. Atomic DOS have been scaled to 
one atom. 



similarly to what happens for the atomic swaps but in 
the case of Cr antisites in NiCrSi we have to reach 20% 
of Cr antisites to loose half-metallicity. Thus for moder- 
ate values of antisites and atomic swaps (around 5-10%) 
both NiCrSi and NiMnSi remain half-metallic although 
the gap shrinks with respect to the perfect alloys. 



B. Defects destroying the half-metallicity 

We will proceed now with the case of defects de- 
stroying the gap and firstly the creation of Ni anti- 
sites at the void site. We have studied two cases 
(i) there is a surplus of Ni atoms at the void site 
[NiCr(Ei_2:Ni2;)Si and NiMn(Ei_2;Ni^)Si compounds 
where E denotes the void site], and (ii) Ni atoms mi- 
grate from the perfect Ni site to the void site conserving 
the number of Ni atoms [(Nii_3;E2:)Cr(Ei_a;Nia;)Si and 
(Nii_j,Ej,)Mn(Ei_2;Nia;)Si compounds]. We present in 
Fig. [T2] the DOS for the Mn-based compounds and for 
a; = 0.1 for both kind of Ni antisites. The results are simi- 
lar for the Cr-based alloys and thus we do not present the 
DOS in this case. We see than in all cases there is a clear 
peak instead of the minority-spin gap. This is not present 
only for the Ni impurity atom but also for the Ni and Mn 
atoms at the perfect sites. In the case of Ni surplus the 
intensity of the peak is larger and is more localized in 
energy. All other details of the DOS are the same with 
the perfect NiMnSi compound and the Ni impurity atom 
has a similar DOS with the Ni atoms at the perfect sites 
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FIG. 13: (Color onlme)FPLO-calculated total and atom re- 
solved DOS for the case of Mn-Ni atomic swaps in NiMnSi 
(left panel) and Cr-Ni atomic swaps in NiCrSi (right panel). 
Atomic DOS have been scaled to one atom. 



since the Void and Ni sites have the same symmetry (four 
Mn and four Si atoms as first neighbors). This peak is 
present also in the case of the NiCrSi alloys and it has the 
same shape and intensity. This gives the hint that this 
peak comes from only one band. To investigate it further 
we show in Fig. [11] also the case of a single Ni impurity 
at a Void site in NiCrSi with the FSKKR impurity code 
(the details are the same with the case of Cr at Si site in 
NiCrSi discussed above). Clearly the peak comes from 
the Cg electrons which form a very narrow band and the 
Fermi level is located exactly at this peak. Contrary to 
the case of a single Cr impurity at a Si site, this peak 
is present for all neighboring atoms and survives practi- 
cally for most of the atoms of the cluster showing that 
the presence of Ni impurity atoms provokes a shift of the 
minority unoccupied antibonding Cg band lower in en- 
ergy and this band has its weight mainly on the Cr(Mn) 
atoms. We have also examined the cases where x= 0.05 
and 0.2 for both NiCrSi and NiMnSi and for both kind 
of Ni antisites. Even for x=0.05 the peak is present but 
now it is more narrow with larger intensity, and as we 
increase the concentration of antisites the peak starts to 
occupy a larger energy range and its intensity becomes 
smaller. Thus we can safely conclude that Ni antisites at 
the vacant site completely kill the spin-polarization. 



TABLE VII: (Color online)FPLO calculated atomic spin mo- 
ments in yLB for the case of Mn-Ni and Cr-Ni atomic swaps in 
NiMnSi and NiCrSi compounds, respectively, within the CPA 
approximation. In the case of the perfect compounds [x = 0) 
we give in parenthesis the spin moment of a single impurity 
atom calculated using the FSKKR method. 



(Nii_^Mn;^)(Mni_:^Ni^)Si 



X 



0.1 
0.2 


Ni 

m 

0.207 
0.180 
0.186 


(0.424) 
0.480 
0.302 


3.005 
2.849 
2.674 


(-1.451) 
-1.496 
-1.358 


m^' 
-0.212 
-0.125 
-0.096 


3.000 
2.509 
1.988 






(Nii- 


-^Cr^)(Cri_^Ni^)Si 




X 


m^' 








m^' 


^total 





0.127 


(0.429) 


2.020 


(-0.870) 


-0.148 


2.000 


0.1 


0.078 


0.409 


1.652 


-0.943 


-0.081 


1.433 


0.2 


0.054 


0.244 


1.196 


-0.656 


-0.052 


0.870 



Finally, we have investigated the case of Ni-Cr and Ni- 
Mn atomic swaps. These defects have the peculiarity that 
Cr(Mn) and Ni atoms move to sites of different symme- 
try with respect to their perfect sites in the NiCrSi and 
NiMnSi alloys. In Fig. [T3]we present for both compounds 
the total DOS for two different concentrations of atomic 
swaps (a;=0.1 and 0.2) and the atom-resolved DOS for 
the first case. The Ni and Cr(Mn) atoms at the perfect 
sites show a DOS similar to the perfect compounds pre- 
sented in Fig. m The impurity Ni atoms at the perfect 
Cr(Mn) sites show a double-peak structure and the Fermi 
level falls within this structure. When the concentration 
is increased to 0.2, as it is obvious from the total DOS, 
the two peaks overlap creating a wide-band. The Ni im- 
purity atoms carry a very large spin moment with respect 
to the spin moment of the Ni atoms at the perfect Ni site 
as can be seen in Table I VIII where we present the spin 
magnetic moments for both Cr- and Mn-alloys and for 
various concentrations of atomic swaps. As we increase 
the concentration of atomic swaps the spin moments of 
both the Ni impurity atoms and Ni atoms at the perfect 
sites decrease. The changes are even more important for 
the Cr(Mn) atoms at the Ni sites. Due to the different 
symmetry these atoms have a spin moment antiparallel 
to the spin moment of the other transition metal atoms 
reducing considerably the total spin moment and leading 
to ferrimagnetic alloys. This behavior is present also in 
the case of Cr and Mn impurity atoms at perfect Co sites 
in Co2CrAl(Si) and Co2MnAl(Si) fuU-Heusler alloys^i22 
where Cr(Mn) impurity atoms have spin moments an- 
tiparallel to the spin magnetic moment of the other tran- 
sition metal atoms. But in these fuU-Heusler alloys the 
half-metallic character is kept contrary to the compounds 
under study here. This is due to the Ni impurity atoms 
and not to the Cr(Mn) impurity atoms which have a very 
small DOS in the minority-spin band which is obviously 
a very weak image of the Ni-impurity DOS. To examine 
this behavior even further we have plotted in Fig. Illl the 
DOS for a single Cr impurity atom at a perfect Ni site 
and for a single Ni impurity atom at a perfect Cr site 
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in NiCrSi. In the case of the Cr impurity atom the gap 
persists and the spin moment is negative (in Table IVIII 
we present in parenthesis the spin moment of the single 
impurity atoms) and close to the values calculated us- 
ing the FPLO-CPA code. The Ni impurity atoms at the 
Cr site show a double peak structure at the Fermi level 
with the Cg states just below the gap and the Fermi level 
being pinned exactly at the minority t2g peak. Thus 
due to symmetry reasons the majority d-states of Ni- 
impurity atoms are completely occupied leading to large 
spin moments and the Fermi level falls within a double 
peak-structure leading to a complete destruction of half- 
metallicity. Thus any kind of Ni defects leads to loss of 
the half-metallic character. 



VII. SUMMARY AND CONCLUSIONS 

We have studied using both the full-potential 
nonorthogonal local-orbital minimum-basis band struc- 
ture scheme (FFLO) and the full-potential screened 
Korringa-Kohn-Rostoker (FSKKR) electronic structure 
methods the electronic, magnetic and gap-related prop- 
erties of the NiYSi compounds and have expanded our 
study also to cover the case of surfaces, interfaces with 
semiconductors and defects. When Y stands for V, the 
ferromagnetism is not very stable due to the weak V- 
V interactions and these compounds are not suitable for 
applications. When Y stands for Cr or Mn, the ferromag- 
netism is extremely stable leading to very high values of 
the Curie temperature which is predicted to be '^700 K 
for NiCrSi and -1100 K for NiMnSi. 

Both NiCrSi and NiMnSi are half-metallic at their 
equilibrium lattice constant with large width of the 
minority-spin gap ('^l eV) and integer values of the total 
spin moment as predicted by the Slater Pauling rule. The 
gap is created due to the creation of bonding and anti- 



bonding d-hybrids in the minority-spin band. The width 
of the gap is marginally affected upon tetragonalization 
even when we expand or contract the lattice by 5%. 

NiCrSi (001) surfaces present a high spin-polarization 
at the Fermi level with respect to the NiMnSi alloy due 
to the large intensity of the Cr majority-spin density of 
states (DOS) at the Fermi level. In the case of inter- 
faces with semiconductors with similar lattice parameter 
(GaP, ZnS and Si), Ni-based contacts show larger spin- 
polarization since Ni atoms at the interface layer have 
a more bulk-like environment and in the case of NiCrSi 
interfaces the high Cr-DOS leads to values of the spin- 
polarization for the Ni-based contacts as high as 90%. 

Finally, we have shown that there are two kind of 
defects. The defects and atomic swaps involving the 
Cr(Mn) and Si atoms lead to a broadening of the bands 
and the gap is slowly shrinking and for a critical value of 
the defects-concentration it disappears. Ni defects on the 
other hand lead to a shift of the energy-localized eg states 
within the gap and these impurity states completely de- 
stroy the half-metallicity. 

For realistic applications it seems that NiCrSi is more 
suitable with respect to both NiMnSi studied here and 
the well-known and widely studied NiMnSb Heuslcr al- 
loy. Its Curie temperature is smaller than the Mn-based 
alloys but it is still high-enough for applications exceed- 
ing considerably the room temperature. Moreover the 
large population of Cr majority-spin states at the Fermi 
level ensures the high-spin polarization at the Fermi level 
which is needed for the injection of current in semicon- 
ductors. Crucial for the operation of devices based on 
Heuslcr alloys is the prevention of the creation of Ni de- 
fects, since they induce impurity states within the gap 
which can couple to interface states and completely de- 
stroy the spin-polarization of the current injected into 
the semiconductor. 
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